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ABSTRACT 

An  experimental  study  was  performed  to  establish  a  technique  for  measuring  the 
volume  reverberation  from  ocean  sediments.   Two  types  of  sediments  (aggregate  and 
fine  sand)  were  used  in  this  study.   The  inhomogeneity  within  the  sediment  caused 
considerable  sample-to-sample  fluctuation  in  the  scattered  waveform.   This  fluctuation 
was  removed  by  spatial  averaging  to  obtain  a  mean  value  over  the  sampling  area.    A 
approximate  model  for  volume  reverberation  from  sediments  was  developed  for  an 
acoustic  pulse  with  an  exponential  decay.   The  results  are  promising.    Combining  the 
model  and  the  experimental  results,  the  volume  scattering  coefficient  obtained  for  the 
aggregate  is  0.0624  ±  0.007  m"3  and  for  fine  sand  is  0.0413  ±  0.007  m"3.    The  latter 
result  is  close  to  the  coefficient  obtained  from  a  cloud-of-small-spheres  model  for 
wavelengths  much  greater  than  the  particle  size. 
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I.  INTRODUCTION 

In  hydrography,  the  character  of  the  bottom  should  be  determined  for  nautical 
charting,  particularly  in  harbors,  anchorage  areas  or  other  areas  where  the  bottom 
characteristic  is  a  significant  factor  to  be  considered  for  navigation  safety. 
Conventional  sampling  methods,  using  clamshell  snappers  (Umbach,  1976).  are  used  to 
define  the  characteristics  of  the  bottom  surface  layer.    It  is  tedious  and  time 
consuming,  especially  in  off-shore  geologic  investigations.    If  some  kind  of  theory  or 
model  existed  to  remotely  classify  the  sediments,  the  task  could  be  carried  out  easier 
and  faster.    Moreover,  this  ability  would  permit  nautical  charts  to  contain  much  more 
information. 

In  the  ocean  bottom,  the  magnitude  and  nature  of  the  sound  scattering  are 
functions  of  both  the  particle  size  of  the  bottom  and  surface  (bottom)  relief  (Mckinney 
and  Anderson,  1964).    The  sound  scho  from  sediments  is  subject  to  three  types  of 
reflecting  or  scattering  (Clark.  Proni.  Seem.  Tsai,  1985)  : 

1)  coherent  reflection  from  the  water-sediment  interface;  the  signal  amplitude 
will  be  approximately  that  given  by  spherical  spreading  of  the  wave  front 
from  the  image  source  modified  by  the  plane  wave  reflection  coefficient 
(Clay,  1966). 

2)  incoherent  or  statistically  variable  scattering  from  the  surface  irregularities 
of  the  sediment,  and 

3)  incoherent  scattering  from  within  the  volume  of  the  sediment  caused  by  the 
acoustically  irregular  matrix  of  the  sediment. 

Items  1  and  2  arise  from  the  surface  of  the  sediments  and  item  3  from  the  sound 
pulse  penetrating  into  the  sediment  volume  and  then  being  scattered  back  to  the 
transducer  by  the  granular  matrix  of  the  sediment.    Since  different  sediments  have 
different  acoustic  matrices,  they  have  different  volume  reverberation  characteristics  and 
there  should  exist  a  mathematical  relationship  between  sediment  properties  and  the 
volume  reverberation. 

The  major  goal  of  this  experiment  was  to  establish  a  fundamental  laboratory 
procedure  for  the  remote  sensing  of  ocean  sediment  volume  reverberation.   This  was  to 
be  accomplished  by  accurately  measuring  the  properties  of  a  pulse  reflected  from  the 
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sediment.    In  application,  since  the  reflected  sound  puise  is  received  by  the  same 
transducer,  i.e.,  "monostatic  sonar',  cur  discussion  is  restricted  to  "back-reflection"  and 
"backscatterine". 
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II.  BACKGROUND 

Bradshaw  (1981)  used  Monterey  #31  fine  sand,  the  same  material  used  in  this 
experiment,  to  study  the  sound  propagation  into  a  fast  bottom  medium.    In  his  study, 
air  bubbles  were  removed  from  within  the  sediment  based  on  heating  the  water  sand 
system.    Bleach  was  added  at  a  ratio  of  one-half  gallon  bleach  to  70  gallons  of  water 
plus  sand  in  order  to  control  biologic  growth.   The  physical  properties,  such  as  density, 
porosity,  speed  of  sound,  reflecton  coefficient  etc.,  were  measured.    Kosnik  (1984)  used 
a  high  pressure  water  jet  to  remove  the  air  bubbles  from  the  water,  sand  mixture. 
However,  no  volume  scattering  measurements  were  made  by  either  investigator. 

Morse  and  Ingard  (1968)  give  a  theoretical  derivation  on  incoherent  scattering. 
Under  the  assumptions  that  :  1)  the  scattering  objects  are  quite  small  compared  with 
wave  length  A,  2)  the  dimensions  of  scattering  region  are  much  larger  than  X,  3)  the 
scatterers  are  all  spheres,  and  4)  the  scatterers  are  populated  sparsely  and  no  mutiple 
scattering  existed,  the  incoherently  scattered  intensity  I-     per  unit  incident  intensity  I-, 
at  angle  9  to  the  incident  wave,  for  frequency  <a,'2n,  is 


Ifc/Ii"  N  V(27rj1/2  K4a6|yK  +yp  cos0|2  (2.1; 


9 


(IK2  a2  sin2(  1/2)0)  exp(-2K2  a2  sin2(  1/2)9) 
8 


where  N  is  the  total  number  of  scatterers  within  unit  volume.  K  is  the  wave  number  for 
the  mixture,  a  is  the  radius  of  scatterers.    Furthermore,  for  monostatic  sonar,  the 
volume  scattering  coefficient  for  unit  volume  can  be  reduced  as 

sy=  X(27i)1'2K6aS|yK-Yp|2  (2-2) 

where 

Kn   "   Ko 


Ko 
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3  P„  "   3  P, 


■P  = 


2  Pn  "   P, 


and  K    is  the  compressibility  of  the  scatterer,  KQ  is  the  compressibility  of  the  fluid,  pn 
is  the  density  of  the  scatterer  and  pQ  is  the  density  of  fluid. 

For  a  small  (ka  <  <   1)  rigid  sphere,  the  backscattering  cross  section  is  given 
(Urick,  1983)  as 

(7  =  2.S  (7C  a2)(K  a)4 

Note  that,  as  a  single  sphere,  it  must  be  a  single  scattering.   The  scattering  coefficient 

for  unit  volume  can  be  obtained  as 

sv  =  N  (J  (2.3) 

Both  theories  are  not  very  appropriate  for  the  real  ocean  sediment  volume 
reverberation.    However,  these  are  the  only  two  theories  available.    The  backscattering 
of  sound  from  the  ocean  bottom  has' been  studied  for  many  years,  but  unfortunately  no 
papers  were  found  to  discuss  the  volume  reverberation  of  ocean  sediments. 
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III.  DESCRIPTION  OF  VOLUME  REVERBERATION 

As  sound  propagates  in  a  fluid,  part  of  the  sound  is  intercepted  and  reradiated  by 
inhomogeneites  in  the  fluid.  The  reradiation  of  sound  is  called  scattering.  The  total  of 
the  scattering  contributions  from  all  the  scatterers  is  called  reverberation  (Urick,  1983). 

Volume  reverberation  is  defined  as  sound  scattered  back  to  the  transducer  by 
scattering  centers  in  the  volume  of  the  sea  (National  Defense  Research  Commitee, 
1969).    Volume  reverberation  may  be  considered  as  a  blending  of  large  number  of 
echoes.    It  has  been  studied  mostly  in  the  Deep  Scattering  layer  ,  a  region  of  high 
biological  activity.    Reverberation  level,  in  decibels,  is  used  to  measure  volume 
reverberation.    Consider  a  directional  projector  in  an  ocean  containing  a  large  number 
of  volume  scatterers  (Figure  3.1).    Let  the  axial  intensity  at  unit  distance  be  I    .   At  a 
distance  r  where  there  is  a  volume  element  dV  of  volume  scatterers,  the  incident 
intensity  is 

The  scattered  intensity  1  m  from  the  volume  dV  is 
dls(l)=  IiSvdV 

where  s,.  is  the  volume  scattering  coefficient.    The  scattered  intensity  back  at  the 
receiver  is 

dls(r)  =  I^dV.'r2 

=  (I0r4)svdV 

and  the  scattered  intensity  from  all  scatterers  contributing  to  the  scattered  intensity  at 
time  t  is 

Is  (r)  =  (I0/r4)  sv  V 

where  V  is  the  total  volume  contributing  for  the  geometry  of  Figure  3.1  calculated  by 
(Urick,  1983) 

V-.Qr2  ex/2 
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Figure  3.1     Portion  of  ocean  scattering. 


IS 


where  H  is  the  solid  angle  of  the  sound  beam.  T  is  the  pulse  length  and  c  is  the  sound 
speed,  so 


Is  (r)  =  (I0  r2)  syn  ct/2 

The  reverberation  level  is 

RL  =  10  LOG  Is(r)/Iref 

=  SL  -  TL  +  Sv  4-  10  LOG  (H  ct/2) 

where  TL  is  the  transmission  loss  for  spherical  spreading  20LOG(r)  and  Sy  is  the 
scattering  strength  for  unit  volume  10LOG(s  ).    In  Deep  Scattering  Layer  studies,  this 
reverberation  level  equation  is  sufficient,  because  the  received  pulse  is  regarded  as 
totally  coming  from  volume  reverberation,  i.e.   negligible  reflection  and  surface 
scattering.    Furthermore,  the  derivation  of  V  is  based  on  a  square  acoustic  pulse.    But 
for  a  real  transducer,  the  square  pulse  is  transformed  into  a  pulse  with  an  exponential 
rise,  a  flattened  top  and  an  exponential  decay.    The  sound  pulse  reflected  from  the 
sediment  and  received  by  the  transducer  combines  coherent  reflected,  incoherent 
scattering  from  the  surface  and  volume  reverberation. 

Let  a  transducer  located  in  the  water  send  out  a  sound  pulse  to  a  smooth 
water  air  interface.   The  received  echo  waveform  is  purely  reflection.    If  the  water 
interface  is  rough,  then  the  received  echo  waveform  includes  both  reflection  and  surface 
scattering.    Since  sound  energy  is  scattered  by  the  rough  surface,  the  received  waveform 
amplitude  is  reduced.    Consider  first-order  scattering  and  suppose  that  the  surface 
roughness  is  the  same  for  every  position.    The  amplitude  of  the  combined  sound  wave 
(reflected  and  scattered)  is  reduced  proportional  to  the  purely  reflected  wave  amplitude. 
In  Figure  3.2(a),  the  solid  line  represents  an  ensemble  of  reflected  waveform,  and  the 
dashed  line  represents  an  ensemble  of  reflected  and  scattered  waveform. 

The  contributions  to  a  pulse  reflected  from  a  sediment  surface  is  shown  on 
Figure  3.2(b).   The  dashed  line  is  the  contribution  by  reflection  and  surface  scattering: 
again,  only  first  order  scattering  is  considered.    The  dotted  line  includes  volume 
reverberation.    The  level  of  the  signal  is  raised  above  that  of  a  dashed  line  because  of 
energy  returned  from  within  the  sediment.    In  particular,  the  tail  decays  more  slowly 
than  for  reflection  from  the  water  air  interface.    This  is  the  significant  feature  of 
sediment  volume  reverberation  that  we  are  after  in  this  paper. 
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WAVEFORM  ENVELOP  FOR  WATER/AIR  SURFACE 


V 


V 


Xs 


WAVEFORM  ENVELOP   FOR  SEDIMENT  SURFACE 


Figure  3.2     Theoretical  echo  waveform  for  water  and  sediment. 
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Two  alternative  assumptions  have  to  be  made  for  reflection  from  the  sediment 
surface  :  either  1)  the  surface  scattering  is  negligible  compared  with  volume 
reverberation,  or  2)  if  it  is  not  negligible,  it  must  have  the  same  amplitude  factor  as 
discussed  in  Figure  3.2(b). 

According  to  these  assumptions,  if  a  sediment  echo  is  subtracted  from  the  water 
echo  it  will  give  a  measure  of  the  volume  reverberation  in  the  sediment  (Figure  3.3). 
The  decaying  tail  of  the  echo  reflected  from  the  smooth  water,  air  interface  is  an 
exponential  decay. 

Vw(t)  =  V0e-* 

where   Vw  (t)  is  the  received  voltage  at  time  t.  5  is  the  decay  constant,  V    is  a  constant 
voltage,  and  t  is  the  time  starting  from  the  beginning  of  decay.    If  we  assume  the 
sediment  echo  also  decays  exponentially,  then  we  will  obtain  a  similar  function 

vs  (t)  =  v'0  e-nt 

where   V   (t)  is  the  received  voltage  from  sediment  at  time  t  and  \\  is  a  decay  constant. 
If  Vw  and  V    are  both  normalized  to  the  same  initial  value  K,  the  difference  becomes 

Vy(t)  =  K  (e-W-e"51)  (3.1) 

If  both  ¥|t  and  5t  are  much  less  than  unity, 

Vy  (t)  =  K  (  11  -  §  )t         t  <  <  1/n  and  1/8 

This  initial  decay  is  linear  and  depends  on  the  difference  S-T|. 
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Figure  3.3     Real  echo  waveform  displayed  on  CRT. 
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IV.  EXPERIMENT  DESIGN 

A.  MATERIAL  AND  TANK  SELECTED 

Fresh  (tap)  water,  #30  Monterey  line  sand  and  Monterey  Aquarium  ^2  for  a 
rough  gravel  sediment  were  the  media  used  in  the  experiment.   The  grain  size  of  #30 
fine  sand  varies  from  0.15  mm  to  0.7  mm  (Bradshaw.  19S1)  with  an  average  size  0.3 
mm.   The  Aquarium  #2  has  a  mean  size  5.3  mm  (Diaz,  1986).   The  sand  water  mixture 
was  repeatedly  stirred  to  remove  air  bubble  trapped  while  transferring  the  sand  to 
different  tanks. 

Initial  experiments  were  conducted  in  the  NTS  anechoic  tank  to  measure  the 
reflection  from  the  water/air  surface,  with  the  transducer  located  at  different  depth. 
The  second  tank,  a  wooden  tank  measuring  SO  cm  x  SO  cm  x  60  cm,  was  constructed  to 
hold  about  30  cm  of  aggregate.    This  tank  was  lowered  and  suspended  85  cm  below  the 
water  air  surface  of  the  anechoic  tank.   A  third  tank,  a  steel  glass  tank,  measuring  70 
cm  x  70  cm  x  60  cm  was  filled  with  24  cm  of  fine  sand  or  aggregate  with  36  cm  of 
water  above. 

B.  ELECTRONIC  EQUIPMENT 

A  schematic  drawing  of  the  equipment  is  shown  in  Figure  4.1.    Except  a  high 
pass  filter,  all  components  were  off-the-shelf.    General  Radio  model  1310  oscillator 
with  a  frequency  of  182.0  kHz  was  fed  simultaneously  into  a  Hewlett-Packard  523 3 L 
frequency  counter,  and  a  General  Radio  Type  396-A  tone  burst  generator  was  used  to 
generate  a  16  cycle  pulse.    Signal  output  from  Hewlett-Packard  IIP  467-A  amplifier 
was  fed  through  a  Datasonic  Transmit  Receive  (T/R)  switch  and  a  small  21  to  2S  kHz 
high  pass  filter,  then  into  the  F-41  transducer.    The  F-41  transducer  was  discussed  in 
detail  by  both  Borchardt  (1985)  and  Diaz  (1986).    Signals  returned  from  the  sediment 
were  amplified  20  dB  by  a  Hewlett-Packard  465-A  pre-amplifier,  then  passed  through  a 
Spencer-Kennedy  Laboratories,  Inc.  model  302  variable  electronic  filter  (set  at  135  kHz 
high  pass)  to  eliminate  low  frequency  noise  before  being  passed  to  the  digital  and 
analog  oscilloscopes. 

The  waveform  displayed  on  3091  Nicolet  Digital  Oscilloscope  was  sampled  and 
sent  as  a  stream  of  characters  to  the  computer.  To  reduce  manual  operations  on  the 
Nicolet  3091  during  the  experiment,  a  IIP-3421 A  Data  Acquisition  and  Control  Unit 
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Figure  4.1     Schematic  of  equipment. 
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was  adopted  to  receive  commands  from  a  HP-S6  computer  and  then  automatically 
operate  the  Nicolet  3091.    However,  it  was  noticed  that  if  the  Nicolet  3091  was  not 
"warmed  up"  enough,  the  DC-offset  jumped  up  and  down,  and  continual  manual 
adjustment  was  needed.    After  2  or  3  hours  of  running,  it  became  stable  and  no  further 
adjustment  was  required. 

C.       VERIFICATION  OF  MEASUREMENTS 

No  matter  how  deep  the  transducer  is  located  below  the  water  air  surface,  as 

long  as  the  transducer  is  level,  the  exponential  decay  of  the  tail  should  be  the  same.    In 

Diaz's  thesis,  the  measurements  were  inconclusive  because  the  transducer  was  fixed  at 

one  position  throughout  the  experiment.    This  was  not  adequate  to  prove  the  system 

functioned  properly.   To  prove  that  the  previous  experiments  were  run  properly  and 

results  were  obtained  correctly,  a  verification  was  attempted  by  measuring  the 

reflection  from  the  water  air  interface  with  the  same  procedure,  but  with  the  transducer 

at  various  depths.    Data  set  20  (Figure  4.2)  for  the  transducer  70  cm  below  the  surface, 

■> 
gives  a  siepe  of  -3.75  x  'J)'~  Np  us.    Data  set  21  (Figure  4.3)  was  obtained  by  raising 

the  transducer  15  cm  and  the  slope  is  -3.77  x  10     Np  us.    Both  results  were  close  to 

the  previous  results  obtained  by  Diaz  (1986).   After  these  two  data  sets  were  collected, 

a  new  transducer  mounting  was  built  and  adopted  to  improve  the  experimental 

precision  and  operation  ability.   The  transducer  was  leveled  by  a  counter-weight  on  the 

new  carriage.   To  ensure  the  system  and  operation  worked  properly,  more  locations 

were  tried  to  obtain  the  decay  constant  of  the  water/air  surface. 

Data  sets  23  and  24  were  obtained  at  the  same  transducer  depth,  but  seperated 

40  cm  horizontally.    Data  sets  24  and  25  were  obtained  at  the  same  horizontal  position 

but  with  50  cm  difference  in  depth.    Data  set  26  was  acquired  by  changing  10  cm  in 

both  the  depth  and  position  from  data  set  25.    Data  sets  26  to  2S  resulted  from 

changing  the  depth  from  70  cm  to  30  cm  below  the  water  air  interface.    The  results  of 

these  experiments  are  listed  in  Table  1  and  show  consistencies  between  each  other. 

Waveform  ensembles  of  the  tail  for  each  data  set  arc  plotted  on  Figure  4.4  to  Figure 

4.9. 

The  average  slope  -3.77  ±  0.077  x  10"2  Np  us  is  smaller  than  -3.S4  ±  0.14  x 
10"^  Np  us  which  was  measured  in  Diaz's  thesis.    Statistically,  these  two  slopes  are  the 
same.   The  smaller  standard  deviation  might  suggest  that  the  new  carriage  has  achieved 
a  better  precision. 
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Figurc  4.2     Decay  envelope  for  reflection  from  water/air  interface, 
transducer  70  cm  below  the  surface  (data  set  20). 
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Figure  4.3     Dccav  envelope  (or  reflection  from  water/air  interface, 
transducer  55  cm  below  the  surface  (data  set  21). 
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TABLE  1 
WATER  AIR  REFLECTION 


DATA 

SET 

DATA20 
DATA2 1 
DATA22 
DATA23 
DATA24 
DATA25 
DATA2  6 
DATA2  7 
DATA28 

MEAN    SLOPE 

c   =   0. 07    x    10 


SLOPE 


(x!0~2) 

-3. 

75 

-3. 

77 

-3. 

86 

~3. 

62 

-3. 

68 

-3 . 

86 

-3. 

80 

-3. 

84 

-3. 

79 

-2 


-3.  77    x    10    *    Np/fis 


CORRELATION 
COEFFICIENT 

-0. 9984 
-0. 9795 
-0. 9979 
-0. 9724 
-0. 9895 
-0. 9956 
-0. 9679 
-0. 9888 
-0. 9810 


-2 


Np/Jis 


After  dealing  with  the  water,  air  surface,  the  transducer  was  turned  over  and 
pointed  downward  toward  the  sediment  which  was  contained  within  the  SO  cm  x  SO  cm 
x  60  cm  wooden  tray.    Different  transducer  positions  were  choosen.    Data  sets  31  to  33 
were  at  the  same  depth,  but  horizontally  3  cm  away  from  each  other.   The  results  show 
the  waveform  were  totally  inconsistant.    Data  sets  33  to  35  were  at  the  same  horizontal 
position  but  different  depths;  results  were  also  inconsistant.   Above  results  are 
presented  on  Table  2  and  plotted  on  Figure  4.10  through  Figure  4.14.    So  far, 
experiments  reveal  the  lack  of  isotropic  homogeniety  of  the  volume  contributes 
reverberation. 

D.       MEASUREMENT  CONSIDERATION 

During  the  course  of  these  experiments  some  other  effects  were  considered  and 
studied  to  ensure  that  the  experiment  was  working  properly. 
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Figure  4.4     Decay  envelope  fur  reflection  from  water/air  interface, 
transducer  was  3d  cm  below  surface  (data  set  23). 
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Figure  4.5     Decay  envelope  for  reflection  from  water/air  interface. 
trnsducer  was  30  cm  below  surface  (data  set  24). 
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Figure  4.6     Decay  envelope  for  reflection  from  water/air  interface, 
transducer  was  Si)  cm  below  surface  (data  set  25). 


31 


UK  (<V0LTAGE"2>)  "1/2) 


m 


r> 


i 
en 

o 

D 


1G30 


1050 


1070  - 


§  1090 
i 

in 
m 
n 


1110 


1130 


1150 


1170 


1190  - 


i 

CO 

o 

o 


I 

o 
o 


p 


o 


o 


en 

o 

r~ 

o 

o 

TO 

"D 

3j 

m 

m 

o 


o 

CO 

O 

03 

CD 

O 

m 

CO 

en 

CO 


Figure  4.7     Decay  envelope  for  reflection  from  water/air  interface, 
transducer  was  70  cm  below  surface  (data  set  26). 
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Figure  4.8     Decav  envelope  for  reflection  from  water/air  interface, 
transducer  was  50  cm  below  surface  (data  set  27). 
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Figure  4.9     Dccav  envelope  for  reflection  from  water/air  interface, 
transducer  was  30  cm  below  surface  (data  set  28). 
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TABLE  2 
SEDIMENT  WATER  REFLECTION  AND  SCATTERING 


DATA 
SET 

DATA3 1 
DATA32 
DATA33 
DATA3  4 
DATA3  5 

MEAN  SLOPE  = 

G   =   0.  85  x  10 


SLOPE 


(x!0~2) 

-2. 

20 

-2. 

98 

-1. 

94 

-3. 

06 

-2. 

90 

-2 


62    x    10" 2    Np/Jis 
Np/Jis 


CORRELATION 
COEFFICIENT 

-0. 9138 
-0. 8908 
-0. 7838 
-0. 8174 
-0. 9844 


1.  Effects  of  Sampling  Rate 

The  sound  signal  being  used  was  fixed  at  1S2  kHz  (5.5  jis  period).   The 
Nicolet  3091  digital  oscilloscope  has  as  its  fastest  sampling  rate  1  jis  which  is 
inadequate  to  give  a  precise  reproduction  of  the  numeric  waveform.    However, 
observations  have  shown  "jitter"  in  the  triggering  of  the  scope,  so  for  consecutive 
pulses  the  scope  would  sample  different  parts  of  the  waveform.   To  obtain  an  accurate 
waveform,  a  large  number  of  waveform  were  sampled  and  averaged.    If  sufficient 
"jitter"  is  present,  the  voltage  measured  in  a  given  bin  should  vary  between  the 
maximum  and  minimum  voltage  in  the  waveform  in  the  vicinity  of  the  corresponding 
time.    Since  the  sine  wave  varies  between  minus  and  plus,  the  voltage  values  were  • 
squared  before  averaging  in  order  to  give  an  accurate  representation  of  the  waveform 
envelope. 

To  determine  the  maximum  number  of  samples  required  for  each  data  set 
(average  waveform),  the  water  air  interface  was  measured  with  different  number  of 
samples.   The  results  showed  that  50  samples  for  each  data  set  were  enough  to  give  a 
true  waveform. 

2.  Effects  of  Bottom  Reflection  Interference 

For  a  16  cycle  signal  with  a  frequency  182  kHz,  the  pulse  length  is  SS  jtis. 
Since  the  sound  pulse  emitted  had  a  decay  time  of  approximately  100  fis,  the  echo  from 
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rigure  4.10    Decay  envelope  for  rcllection  from  aegregate/water 

interface,  transducer  was  30  cm  above  surface  (data  set  31). 
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Figure  4.1 1     Decay  envelope  for  reflection  from  aggregate/water 

interface,  transducer  was  30  cm  above  surface  (data  set  32). 
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Figure  4.12     Decay  envelope  for  reflection  from  aggregate/water 
interface,  transducer  was  30  cm  above  surface  (clata  set  33). 
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Figure  4.13     Dccav  envelope  for  reflection  from  aggregate/water 
Interface,  transducer  was  50  cm  above  surface  (data  set  34). 
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I-igure  AAA    Decay  envelope  lor  reflection  from  aggregate/water 
Interface,  transducer  was  70  cm  above  surface  (data  set  35). 
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the  smooth  water  air  interface  should  also  have  the  same  decay  time.    Since  the 
purpose  of  this  project  was  to  compare  the  decay  portion  of  echos  reflected  from  water 
and  sediments,  the  waveform  was  sampled  for  a  100  us  interval  from  the  begining  of 
decay.   The  sound  speed  in  the  aggregate  was  measured  to  he  1555  m  s.    For  the  fine 
sand.  Bradshaw  (1981)  obtained  the  speed  of  sound  to  be  1610  m  s.    Combining  these 
sound  speeds  and  the  time  interval  1SS  us,  the  depth  of  sediment  that  contrubutes  to 
the  main  echo  is  14.6  cm  for  aggregate  and  15.1  cm  for  fine  sand.    Interference  caused 
by  reflection  from  the  bottom  of  the  tank,  will  not  affect  the  experiment  as  long  as  the 
sand  is  thicker  than  15.1  cm  and  the  aggregate  is  thicker  than  14.6  cm.    Throughout 
our  experiments,  we  kept  both  sediments  more  than  24  cm  deep  in  the  glass  tank  to 
avoid  interference. 

3.  Effects  of  Near  Field  Interference 

In  the  near  field  of  a  transmitter,  the  axial  pressure  exhibits  strong  interference 
effects  (Kinsler.  Frey.  Coppens.  Sanders,  19S2).    Thus,  the  sound  field  is  irregular  and 
does  not  fall  off  smoothly  with  distance  as  it  does  in  the  far  field.    To  avoid 
measurements  in  the  near  field,  an  experiment  was  run  to  find  the  transition  distance 
which  separates  the  near  field  and  far  field.   This  distance  is  20  cm  from  the  transducer 
center.   No  problems  should  occur  as  long  as  the  target  is  more  than  20  cm  from  the 
transducer.    In  the  glass  tank,  we  placed  the  transducer  30  cm  above  sediment.   The 
water 'air  interface  measurement  also  gave  a  check  that  30  cm  distance  is  satisfied. 
Since  a  pulse  moving  in  water  at  14S0  m  s,  takes  405  us  to  travel  the  30  cm  to  the 
sediment  interface  and  back,  the  30  cm  distance  is  long  enough  so  that  the  transducer, 
which  is  on  about  1SS  us,  stops  ringing  before  the  first  echo  from  the  sediment  arrives. 

4.  Effects  of  Water/ Air  Interface  Interference 

Ideally,  the  transducer  should  be  far  enough  away  from  all  unwanted  surfaces 
to  avoid  interference.   As  the  transducer  faces  down  and  projects  a  major  part  of  its 
energy  to  the  sediment  surface,  some  energy  is  transmitted  to  the  water  surface  and 
reflected  back  toward  the  sediment.    If  the  transducer  is  close  to  the  water/air  surface, 
the  reflected  energy  will  interfere  with  the  transmitted  energy.    Restricted  by  the  size  of 
the  glass  tank,  we  had  to  locate  the  transducer  close  to  the  water  surface.    To  test  how 
important  this  interference  could  be,  the  water  surface  was  agitated  by  fingers.   The 
waveform  displayed  on  both  analog  and  digital  oscilloscope  was  unchanged  during  the 
ruffling.    It  was  concluded  that  this  source  of  interference  is  negligible  for  this 
experiment. 
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5.  Effects  of  Side  Lobe 

The  half  beamwidth  (to  the  first  null)  for  the  F-41  transducer  was  measured  to 
be  10.0°.    For  a  30  cm  distance  from  sediment  to  the  transducer,  the  sediment  surface 
ensonified  by  the  major  lobe  is  a  circle  area  with  radius  5.3  cm.    If  the  distance  is  SO 
cm.  the  radius  will  be  14.1  cm.    For  both  the  glass  tank  and  wood  tank,  if  the 
transducer  is  located  above  the  center  of  the  tank,  the  major  lobe  will  not  produce 
interference  caused  by  the  sides  of  the  tank. 

The  glass  tank  was  too  small  to  avoid  the  interference,  if  any,  caused  by  the 
side  lobe.    To  test  how  much  the  interference  was,  sediments  outside  the  center  30  cm  x 
30  cm  square  area  were  disturbed  and  left  in  a  very  irregular  surface.    Waveforms 
returned  from  sediment  were  measured  before  and  after  the  disturbance.    Data  set  51 
(Figure  4.15)  shows  the  waveform  measured  before  the  sediment  was  disturbed  and 
data  set  52  (Figure  4.16)  shows  the  waveform  after  disturbance.    By  comparing  the  two 
waveforms,  it  is  found  that  the  F-41  trandsucer  gave  negligible  side  lobe  interference  in 
this  experiment. 

6.  Effects  of  Sediment  Inhomogeneity 

Before  the  test  findings  on  side  lobe  effects,  we  were  bothered  by  the  drastic 
changes  on  the  waveform  whenever  the  transducer  was  moved  parallel  to  the  smoothed 
aggregate  surface.    It  was  assumed  that  the  waveform  fluctuation  was  caused  by  the 
inhomogeneities  of  the  aggregate.   To  prove  this,  we  placed  the  transducer  at  a  fixed 
position  in  the  glass  tank  and  repeatedly  ensonified  the  aggregate.   After  each  data  set 
was  acquired,  the  aggregate  was  disturbed  and  smoothed  again.    Ten  sets  of  waveform 
data  were  obtained  and  each  of  them  differed  from  the  others. 

It  was  thought  the  fine  sand  would  be  more  homogeneous  than  aggregate  and 
it  might  produce  more  consistant  waveforms  as  the  transducer  was  moved  horizontally. 
Measurements  made  with  the  fine  sand,  however,  showed  that  the  waveform  still 
fluctuated.   After  each  data  set  was  obtained  from  the  smoothed  fine  sand,  the 
transducer  was  moved  horizontally  to  different  position.   All  the  transducer  positions 
were  at  least  15  cm  from  the  wall  to  avoid  interference  caused  by  reflection  from  the 
sides  of  tank.   Also,  the  positions  were  at  least  a  length  of  tranducer  diameter  from 
each  other  to  avoid  a  waveform  replication.   A  metal  plate  was  set  on  top  of  the  sand 
and  the  transducer  was  moved  horizontally.   As  expected,  the  waveform  on  the 
oscilloscope  was  unchanged  but  totally  different  from  sand. 
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Figure  4.1.5     Decay  envelope  for  reflection  from  aggregate/ water 
interlace  in  the  glass  tank  (data  det  5  r)T 
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E.        PROGRAMMING 

An  HP-S6  desktop  computer  was  adopted  to  collect  and  analyze  the  waveform 
data.    The  computer  program  was  written  to  he  interactive.    This  enables  1)  the 
collection  of  data  sets  with  different  delay  time  intervals  and  2)  the  collecting  and 
processing  of  numerous  data  sets.    The  interactive  programs  allowed  essential 
information  to  be  entered  via  the  key  board,  eliminating  the  need  to  modify  the 
program  whenever  data  sets  were  changed.    The  interactive  prompts  were  useful  for 
people  unfamiliar  with  Basic  programming. 

1.  Program  "THESIS3" 

This  program  was  designed  as  the  initial  step  in  the  data  collection  processing. 
Before  the  selected  waveform  data  displayed  on  the  3091  Nicolet  Digital  Oscilloscope 
were  output  to  the  computer,  it  was  switched  into  "freeze"  or  "store"  mode  from  "live" 
mode.    Through  the  interface  RS-232.  the  computer  received  the  stored  waveform  data 
and  then  processed  them.    After  processing  the  waveform  data  and  before  receiving  the 
next  sampling  data,  the  Xicolet  3091  had  to  be  switched  into  "live"  mode  again.    An 
HP-342IA  Data  acquisition  and  Control  Unit  was  interfaced  between  the  Nicolet  3091 
and  HP-S6.   The  HP-3421A  reduced  the  amount  of  manual  operations  in  the  data 
collection  process.   A  five  volts  source  were  input  to  the  HP-3421A.    Commands  form 
the  HP-S6  and  5  volts  signal  were  used  to  control,  by  closing  and  opening  the  channel 
in  the  HP-3421A,  the  "mode  switching"  of  the  oscilloscope.   To  avoid  errors,  the 
oscilloscope  should  be  warmed  up  enough  while  entering  the  inital  information.   Once 
this  initial  information  was  entered  and  the  oscilloscope  was  warmed  up  enough;  the 
whole  system  runs  automatically. 

Waveform  data  output  from  the  oscilloscope  were  a  stream  of  characters.    To 
obtain  the  voltage  and  time  value,  all  these  characters  needed  to  be  normalized  and 
converted  by  the  computer.    A  program,  listed  in  the  Nicolet  3091  operation  manual. 
was  used  for  data  value  normalization.   A  data  set,  the  waveform  of  the  decay  portion. 
was  obtained  by  squaring  and  averaging  50  samples  of  waveform  data  for  each 
transducer  position.   The  running  average  formula  used  in  the  program  was 

<  V2>n  =  ((n-l)<  V2  >n.j  +Vn2)n 

where    <  Vz  >     is  n  samples  averaged  squared-voltage,  <  Vz  >n.j  is  n-1  samples 
averaged  squared-voltage  and  V  "  is  squared  voltage  for  the  nth  sample.   Note  that 
the  average  window  for  the  waveform  is  100  fis  from  the  beginning  of  the  decay.    The 
above  formula  was  processed  for  each  corresponding  time  bin. 
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During  each  data  set  processing,  some  error  messages,  such  as  STRING  OVF 
(which  stands  for  the  string  variable  overflow),  were  displayed  on  the  HP-S6  screen. 
This  error  messaged  could  be  fatal  to  the  data  set  if  the  data  has  not  been  stored  on 
the  disk.   To  avoid  this  problem,  the  program  was  designed  so  that  after  every  10 
samples  processed  the  updated  data  were  stored  on  the  disk  over  the  previous  updated 
data.   The  program  eliminated  the  problem  of  the  disk  being  full  before  the  final  data 
set  was  stored  on  the  disk.    This  prevention  could  increase  the  processing  time,  but 
compared  to  the  time  wasted  from  the  above  problems,  it  is  worthwhile.    The  final 
result  for  this  program  was  3  data  files  :  one  stored  the  squared  voltage,  one  for  time 
data  and  one  for  the  number  count  of  the  samples. 

2.  Program  "MAXIMUM" 

In  order  to  obtained  the  envelope  of  the  waveform,  this  program  was  designed 
to  pick  up  the  local  maximum  of  the  voltage  data.   The  voltages,  stored  in  the  data  file 
from  "THESIS3",  were  read  and  compared  with  the  neighboring  voltage  data  to  decide 
the  local  maximum.    The  determined  local  maximums  and  corresponding  time  bins 
were  stored  in  seperated  files  and  printed  out.    The  data  Hies  were  differentiated  by 
number  according  to  the  order  they  were  acquired. 

3.  Program  "PLOT3" 

After  the  local  maximums  were  determined,  the  decay  portion  of  the  waveform 
can  be  plotted  by  program  "PLOT3".   The  experiments,  for  this  project,  took  place 
mainly  in  the  glass  tank.   A  decay  envelope  measured  in  the  glass  tank  has  an  interval 
between  520  and  619  ).is  starting  from  transmitting  time.    Some  exceptions  existed  and 
gave  different  time  frames. 

To  determine  the  decay  constant,  the  voltage  data  were  converted 
logarithmically  and  plotted  vs.  the  time 

V  =  V0  e"at 

LN  (  V  )  =  -at  +  LN  (  VQ  ) 

If  the  cursor  on  the  oscilloscope  was  not  centered,  the  resulting  envelope  showed  a 
"zig-zag"  instead  of  a  straight  line  for  a  perfect  exponential  decay. 

The  least-square  method  was  adopted  to  calculate  the  best  fit  of  straight  line 

y  =  ax  +  b 
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where  a,  which  is  similar  to  the  decay  constant,  is  the  slope  and  b  is  the  intercept  on 
the  y-axis.   According  to  Beers  (1953), 

kl(xnyn)-lxnyyn 


kVx2-(V 


^  -n       vL  *n  ) 


—  An    —  -  n     —  An~  ^An-  n' 


b  =■ 


kV  xn2  -  (V  x 


where  k  is  the  total  number  of  pairs  of(x.y)  values,  n  is  the  number  between  1  and  k, 

y  x_,  is  the  sum  ot  the  total  x  value.  Y  y„  is  the  sum  of  the  total  y  value.  V  x„    is  the 
<—    n  •  *—  *  n  v  — -    n 

sum  of  the  squared  x  value,  y  yRz  is  the  sum  of  the  squared  y  value  and  ^  (xnyn)  1S 
the  sum  of  the  multiples  of  x    and  y 

Note  that  Y  x  ~  and  y  (x„)    are  not  the  same.   Y  x  ■    implies  that  each 
value  of  x  is  squared  and  then  a  sum  is  made  of  these  squares,  while  Y  (xn)    implies 
that  all  of  the  values  of  x  added  together  and  then  this  sum  is  squared.   There  is  an 
analogous  distinction  between  Y  (x„v,.>  and  y  x„y  v„ .    In  our  case,  the  x-coordinate 
is  the  time  value  and  the  y-coordinate  is  the  natural-log  of  the  voltage.   To  check  how 
close  the  observed  and  least-square  estimated  y-value,  compare,  a  correlation 
coefficient  was  calculated  (Chatfield,  19S4) 

(  xn  -  x  )(  yn  -  y  ) 


(  I  (  xn  -  x  )2X  (  yn  -  y  )2  )12 

where  r  is  the  correlation  coefficient,  x  is  the  mean  value  of  x  and  y  is  the  mean  value 
of  y.    Slope  and  correlation  coefficient  were  both  printed  with  the  graph. 
4.  Program  "AVERAGE" 

The  received  echo  waveform  varies  from  position  to  position  because  of  the 
inhomogeneous  sediment.   To  find  the  decay  constant,  we  have  to  average  out  the 
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inhomogeneities  to  obtain  a  mean  value  of  the  waveform.   This  program  performs  a 
running  average  on  the  original  data  sets,  producing  the  final  averaged  waveform  of  the 

t. 

The  above  four  programs  are  given  in  Appendix  A. 


48 


V.  RESULTS 

The  data  sets  for  reflection  from  aggregate  water  interface  are  data  42  to  51. 
Figure  C.l  to  Figure  C.9,  and  Figure  4.15.    Table  3  shows  the  slope  obtained  for  each 
data  set  varied  between  data  sets.    To  get  a  mean  value  for  the  slope,  the  running 
average  method  was  used  to  process  the  da:.;  sets.    The  running  average  method  also 
shows  how  fast  the  average  processing  became  stable.    Data  sets  DAT2  to  DAT  10. 
Figure  CIO  to  Figure  CIS,  are  the  results  for  each  average  process.   The  slopes  are 
listed  on  Table  4.    It  is  evident  that  the  slope  becomes  stable  after  8  sets  of  data  are 
averaged.   The  final  averged  slope  is  2.13  x  10     Np  jis  and  is  quite  different  from  the 
value  reflected  from  the  water  air  interlace.    The  dilTerence  between  the  two  slopes 
resulted  from  the  effects  cf  volume  reverberation  in  the  aggregate,  and  the  scattering  is 
significant. 

TABLE  3 
ORIGINAL  AGGREGATE  REFLECTION 


DATA 
SET 

DATA42 
DATA43 
DATA44 
DATA45 
DATA46 
DATA47 
DATA48 
DATA49 
DATA 50 
DATA51 


SLOPE 
(xlO"2) 


-0.  94 
-1.  35 
-2.  68 
-2.  64 
-0.  12 
-1.  81 
-0.  95 
-3.  77 
-2.  47 
-0.  16 


CORRELATION 
COEFFICIENT 

-0. 5273 
-0. 5874 
-0. 7848 
-0. 8509 
-0. 1162 
-0. 3971 
-0. 3990 
-0. 8654 
-0. 9756 
-0. 0584 


The  data  sets  for  reflection  from  Fine-sand  water  interface  are  data  55  to  64. 
Figure  D.l  to  Figure  D.10.   The  slope  of  each  data  set  is  listed  on  Table  5.    Evidently. 
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TABLE  4 
AVERAGED  AGGREGATE  REFLECTION 


DATA 

SET 

DAT2 
DAT3 
DAT4 
DAT5 
DAT  6 
DAT7 
DATS 
DAT9 
DAT  10 


SLOPE 
(xlO"2' 

-1.  19 

-1.  97 
-2.  32 
-2.  18 
-2.  10 
-1.  98 
-2.  16 
-2.  20 
-2.  13 


CORRELATION 

COEFFICIENT 

-0. 6C51 
-0. 9720 
-0. 9175 
-0.  9343 
-0. 9914 
-0. 9952 
-0. 9932 
-0. 9943 
-0. 9958 


the  waveform  changes  very  much  for  different  data  sets.   The  same  running  average 
method  was  used  to  remove  the  inhomogenities.   Data  sets  DAT  12  to  DAT20  were 
produced  from  the  average  process.   They  are  listed  on  Table  6  and  plotted  from 
Figure  D.l  1  to  Figure  D.19.    The  slope  became  stable  after  6  data  sets  were  averaged. 
This  might  suggest  that  the  fine  sand  has  a  greater  homogeneity.    The  averaged  slope 
for  the  fine  sand  is  2.37  x  10     Np,  ]is  and  is  greater  than  the  aggregate.    However,  the 
slope  difference  between  them  is  not  as  distinguished  when  compared  to  the  slope 
difference  between  the  sediment  and  the  water  air  surface. 

According  to  Equation  3.1,  the  decay  portion  volume  scattering,  represented  in 
voltage,  for  the  aggregate  is 

Va=  K(e-21300t   -   e-37700t) 

and  for  the  fine  sand  is 


V    =  K  {  e-2370ut   -   e" 


37700t 


The  above  relations  only  existed  without  surface  scattering  or  the  surface  scattering 
was  negligible. 
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TABLE  5 
ORIGINAL  FIXE  SAND  REFLECTION 


DATA 
SET 

DATA55 
DATA56 
DATA57 
DATA58 
DATA59 
DATA 60 
DATA61 
DATA62 
DATA63 
DATA64 


SLOPE 
(xlO~2) 

-1.  51 
-3.  03 
-2.  50 
-0.  75 
-1.  00 
-3.  41 
-3.  33 
-1.  42 
-2.  41 
-1.  03 


corrslatic:; 

coefficient 

-0. 8173 
-0. 9581 
-0. 8112 
-0. 7290 
-0. 6776 
-0. 9957 
-0. 9050 
-0. 4922 
-0. 9389 
-0. 5784 


A  simplistic  model  for  sediment  volume  reverberation  was  developed  and  is 
described  in  Appendix  A.    Combining  the  experimental  results  and  the  newly  developed 
model,  the  volume  scattering  coefficient  obtained  was  0.0413  ±  0.007  m     for  fine  sand 
and  0.0624  ±  0.007  m     for  aggregate.    These  value  are  reasonable  compared  to  the 
difference  between  the  decay  constant  obtained  from  the  experiment.    They  are 
produced  as  a  best  fit  using  the  least  square  method.   The  correlation  coefficients, 
approximately  0.8S,  indicate  they  are  quite  a  good  fit.   As  a  first  model,  the  result  is 
promising.    However,  the  discrepancy  between  the  experiment  and  model  might  suggest 
the  model  is  still  very  crude.   Also  mentioned  in  the  Appendix  A  is  that  this  is  only  a 
first  step  to  try  to  obtain  a  theory  to  describe  the  volume  scattering   effect  of  the 
sediment;  it  need  to  be  improved  by  further  effort. 

The  ka  value  for  fine  sand  is  0.12,  which  is  quite  small  compared  to  1.    From 
Equation  2.2.  cloud-of-small-spheres  model,  the  volume  scattering  coefficient  for  the 
fine  sand  is  0.0082  m     and  is  not  too  far  from  0.0413  m    .   As  mentioned  in  Chapter 
2,  Equation  2.2  was  derived  under  various  assumptions  and  does  not  apply  to  real 
sediment.    From  Equation  2.3,  single-small-spherc  model,  the  volume  scattering 
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TABLE  6 
AVERAGED  FIXE  SAND  REFLECTION 


DATA 
SET 

DAT  12 
DAT  13 
DAT  14 
DAT  15 
DAT  16 
DAT  17 
DAT  18 
DAT  19 
DAT20 


SLOPE 
( xlO"2) 

-1.  84 

-2.  09 
-1.  46 

-1.  47 
-2.  36 
-2.  39 
-2.  39 
-2.  38 
-2.  37 


CORRELATION 
COEFFICIENT 

-0. 9265 
-0. 9051 
-0. 9309 
-0. 9223 
-0. 9912 
-0. 9904 
-0. 9902 
-0. 9903 
-0. 9916 


coefficient  is  0.34  m     which  is  an  order  of  magnitude  different  from  the  the  measured 
value  0.0413  m    .   These  values  suggest  that  the  cloud-of-small-spheres  model  is  closer 
to  the  newly  developed  model  than  the  single-small-sphere  model. 

The  scattering  coefficient  obtained  from  Equation  2.2  is  9456  m     and  from 
Equation  2.3  is  1134  m    .    Obviously  the  two  values  deviate  considerably  from  the  true 
value.   The  ka  value  is  equal  to  2.1,  which  is  outside  the  Rayleigh  scattering  limit. 
Thus  Equations  2.2  and  2.3  can  not  be  applied  to  the  aggregate. 
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VI.  CONCLUSION  AND  RECOMMENDATION 

The  following  conclusions  are  possible: 

1.  The  inhomogenkies  within  the  sediment  caused  the  reflected  waveform  to 
change  considerably  for  various  transducer  positions.    The  changes  can  be 
removed  by  averaging  over  the  inhomogeneities  and  the  averaging  yields  a 
mean  decay  constant. 

2.  The  glass  tank  works  very  well,  and  therefore  a  bigger  tank  to  contain  the 
sediment  is  not  necessary. 

3.  The  laboratory  procedure  to  measure  the  volume  reverberation,  using 
off-the-shelf  equipment  and  a  high  pass  filter  in  a  water-over-sediment 
system,  is  workable. 

4.  The  basic  program,  written  for  a  HP-86  desktop  computer,  is  sufficient  for 
:h:s  experiment. 

5.  Although  Equations  2.2  and  2.3  can  not  be  applied  to  aggregate  since  it 
has  a  ka  greater  than  1,  the  aggregate  does  have  a  higher  volume 
reverberation  than  fine  sand.   The  aggregate  has  a  larger  size  than  the 
fine  sand;  this  might  suggest  that  the  larger  size  sediment  causes  more 
volume  reverberation. 

6.  From  the  scattering  coefficient  0.0413  ±  0.007  m     for  the  fine  sand  and 
0.0624  ±  0.007  m     for  the  aggregate,  the  scattered  energy  received  by 
transducer  is  quite  small  compared  to  the  incident  energy. 

7.  For  fine  sand  sediment,  the  cloud-of-small-spheres  model  works  better 
than  the  single-small- sphere  model. 

8.  The  simplistic  model,  derived  from  the  true  acoustical  pulse  and  described 
in  Appendix  A,  gives  a  crude  formulated  volume  scattering  relationship. 
Combining  with  the  experiment  result,  the  sediment's  volume  scattering 
coefficient  can  be  estimated. 

Further  study  of  the  volume  reverberation  from  ocean  sediment  is  recommended. 
Future  experimentation  should  include  further  variation  of  the  sediment  size  and  of  the 
acoustic  frequency  to  find  a  relationship  combining  the  ka  value  and  volume 
reverberation.   The  developed  model  should  be  improved  to  obtain  a  better  result 
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coincident  with  the  experimental  result.   The  absorption  coefficient  of  different 
sediment  with  the  frequency  used  should  be  measured  to  compare  the  value  obtained 
from  Urick's  formula  described  in  appendix  A.    Equations  2.2  and  2.3  should  be 
improved  to  obtain  the  scattering  coefficient,  which  can  be  used  to  check,  the  above 
model.    Continued  use  of  the  experimental  procedure,  equipment  and  computer 
programs,  written  for  this  project  on  the  HP-S6,  is  highly  recommended. 
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APPENDIX  A 
DERIVATION  OF  SCATTERING  FROM  THE  SEDIMENT 

The  generated  electrical  pulse  can  be  described  by  the  step  unit  function,  H,  as 
follows: 

I(t)  =  IQ(  H(t)-H(t-t))  (A.l) 

where  I(t)  is  the  intensity  at  time  t,  IQ  is  the  intensity  amplitude,  t  is  the  pulse  length, 

0  t  <  0 
H(t)  = 

1  t  >  0 

and 

0  t-T  <  0 
H(t-t)  = 

1  t-T  >  0 

Since  the  transducer  is  a  resonant  system  with  damping,  the  square  electrical  pulse  was 
transformed  into  an  acoustic  pulse  with  an  exponential  rise,  an  approximately  flattened 
top  and  exponential  decay.    The  shape  of  the  acoustic  waveform  can  be  described. 
(Sanders  and  Coppens,  19S6),  as 

f(t)  =  (1  -  e"at)  H(t)  +  (1  -  e"aT  )  Il(t-i)  (A.2) 

where  t  is  the  delay  time  and  a  is  the  decay  constant. 

Let  a  transducer,  located  in  the  water,  send  out  a  sound  pulse  to  a  smoothed 
sediment  surface.    To  assure  the  surface  scattering  is  negligible  and  the  volume 
scattering  is  only  single  scattering,  reflected  pulse  received,  by  the  same  transducer, 
should  have  an  intensity  (Sanders  and  Coppens,  1986) 

I  =.(J  Tj2  sv  eft  I0  f(t-2^/c)d^/(r+^)2  +  R{  I(t)  r2  )  /  r2  (A.3) 
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Where  T:  is  the  intensity  transmission  coefficient.  sy  is  the  volume  scattering 
coefficient,  P  is  the  absorption  coefficient  in  the  sediment,  £  is  the  affected  sediment 
depth.  R:  is  the  intensity  reflection  coefficient  and  r  is  the  distance  between  the 
transducer  and  sediment.    From  above  equation,  the  intensity  contributed  by  scattering 
is 


-.2,    ,-P4 


->    -> 


Is  =  jf-  sy  e"^  I0  f(t-2Vodc  (T+tf  r  (A.4) 

and  the  intensity  contributed  by  the  reflection  is 

Ir  =  Rj  I0  f(t)  r4  (AS) 

If  r  >  >  %.  then 

Is  =  (l-2  sy  I0  r4  )jyf^  f\t-2>  c)d£  (A.6) 

where 

f(t-2^c)  =  (1  -  q<1-2^c)  )  H(t-2£  c)  -  (1  -  e-a^T"2^/c)  )  H(t-T-2£/c)  (A.7) 

During  the  integration,  the  lower  boundary  of  c,  is  zero;  the  upper  boundary  is  ct  2  for 
H(t-2£,/c)  and  c(t-T)  2  for  H(t-T-2£,,'c).    Note  that  only  the  decay  portion,  t  >  r,  was 
interested.   After  the  integration, 

Is  =  A(  -(e-cPt/2  -  l)/p  -  (e-cPt/2"  e-at)(-p  +  2a  c)  (A. 8) 

+  (e-^-1)  -  i)/p  +  (C-cP(t-T);2  .  e-a(t-x)  )/(.p  +  2a/c)) 

where 

A  =  Tx2  sy  IQ  /  r4  (A.9) 

From  the  experimental  formula  P  =  0.25f  (Urick,  1979)  and  f  is  182  kHz,  the  absorption 
coefficient 
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P  =  45.5  dB,m  =  5.2  Np/m, 

cp/2  is  4030  Np  s  for  aggregate  and  is  41S6  Np/s  for  the  fine  sand.  The  intensity  decay 
constant  a  is  twice  the  pressure  decay  constant  5.  From  the  experiment,  6  is  37700 
Np/s  and  it  gives  an  a  equal  to  75400  Np/s;  2a  c  is  97.2  Np  m  for  the  aggregate  and  is 
93.6  Np  m  for  the  fine  sand.  After  plugging  in  the  above  values  for  the  aggregate, 
Equation  A.S  becomes 

Is  =  A(e-4030t    11.6)  -  A(S.3  e-75400t)  (A.10) 

For  the  fine  sand,  Equation  A. 7  becomes 

Is  =  A(e'41S6t    11.0)  -  A(S.6  e'75400t)  (A.  11) 

Rj  =  R2 
Ri  +  Ti  =  1 

The  pressure  reflection  coefficient  R  is  0.356  for  the  aggregate  (Diaz,  19S6)  and  0.36 
for  the  fine  sand  (Bradshaw,  1981).    For  the  aggregate, 

Rj  =  0.1270 

Tj  =  O.S730 

For  the  fine  sand, 

Rj  =  0.1296 
T|  =  0.S704 

The  pulse  length  is  88  fis  and  from  equation  A. 5, 

Ir  =  ^R^e-^-De-754001^4  (A.12) 

Using  the  above  values,  the  reflected  intensity  for  the  aggregate  becomes 

Ir=  96.6e-75400tlo/r4  (A.  13) 

and  for  the  fine  sand 

Ir=  98.5  e-75400t  I0 /r4  (A.  14) 
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The  intensity  contributed  by  the  scattering  for  the  aggregate  is 

Is  =  (sv  IQ   r4  )(  e-4030t,  15.2   -   6.3  e"754wt  )  (A.  15) 

and  for  the  fine  sand. 

Is  =  (sv  I0   r4  )(  e"41S6t  14.5   -   6.5  e"75400*  )  (A. 16) 

Combining  the  above  equations,  the  reflected  intensity  from  the  aggregate  received  by 
the  transducer  is: 

I  =  (I0r4)B  (A.17) 

where 

B  =  (  sv  e"4030t  15.2  -   6.3  sy  e'754001  +  96.6  e-7^400t  )  (A. IS) 

For  the  fine  sand, 

I  =  (I0'r4)D  (A.19) 

where 

D  =  (  sy  e-41S6t,  14.5   -   6.5  sy  e"754001  +  9S.5  e'754001  )  (A.20) 

1  '7 
By  inputting  a  series  of  scattering  coefficients  and  computing  the  Log  (B)       and 

Log  (D)  '"",  the  theoretical  envelope  was  determined.    A  least-square-fit  slope  was 

obtained  for  each  envelope.   The  scattering  coefficient  for  each  sediment  can  be 

estimated  by  the  value  of  s     for  which  the  slope  of  the  theoretical  envelope  is  equal  to 

the  slope  of  experiment  data.    For  the  fine  sand,  measured  the  decay  constant  23700 

Np/s  gave  a  scattering  coefficient  of  0.0413  m    .    For  the  aggregate,  the  scattering 

coefficient  0.0624  m     was  obtained  from  the  measured  decay  constant  21300  Np  s. 

Figure  A.l  and  Figure  A. 2  show  the  envelopes  determined  theoretically  from  the 

above  values  of  sy  for  the  aggregate  and  fine  sand  respectively.   The  theoretical 
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envelope  is  plotted  by  normalizing  the  left  end  of  the  best  fit  slope  line  to  the 
value  obtained  from  the  experiment  and  then  apply    g       s  normalization  to  t    . 
theoretical  envelope. 

Results  from  the  fine  sand  were  selected  tc  find  the  uncertainty  of  scattering 

ft 

coefficient  obtained  from  the  above  method.    Bye::-;.--   ring   l  series  of  envelope  plots 
foi  different  scattering  cceiTicie:::.  the  uncertainty  ±  7  was  determined     Figure  A. 3 

and  A.4  show  the  envelopes  for  scattering  coefficients  —  7  different  from  0.  413 

m~J.    The  slope  changes  as  the  scattering  coefficient  changes.    In  the  beginning  of 
decay,  the  dominant  term  in  equation  A.20  is  the  reflection  term  which  is  independent 
on  sv,  and  the  scattering  term  which  depends  en  s,,  is  negligible.   The  theoretical  value 
will  not  change  much  in  the  beginning  of  decay  as  the  scattering  coefficient  char.- 
The  normalization,  described  previously,  still  gave  a  good  approximation  to  the  plot  of 
the  theoretical  envelope  vs.  the  experimental  envelope  ford..;"::;:.:  scattering 
coefficients. 

Note  that  the  m  )del  derived  and  described  in  .his  appendix  is  crude.   A  latei 
improvement  on  this  model  is  suggested. 
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Figure  A.l     For  the  aggregate,  the  decay  envelope  obtained  from  the  model 
vs.  the  decay  envelope  obtained  from  the  experiment. 
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Figure  A. 2     For  the  fine  sand,  the  decay  envelope  obtained  from  the  model 
vs.  the  decay  envelope  obtained  from  the  experiment. 
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Figure  A. 3     For  the  fine  sand,  the  decay  envelope  obtained  from  the  model 
vs.  the  decay  envelope  obtained  from  the  experiment. 
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Figure  A. 4     For  the  fine  sand,  the  decay  envelope  obtained  from  the  model 
vs.  the  decay  envelope  obtained  from  the  experiment. 
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APPENDIX  B 
COMPUTER  PROGRAM 


•  000 
1 0  1 0 
1020 
1 030 
1040 
1050 
1  060 
1070 
1 080 
1090 
1  1 00 
1110 
1  120 
1  1  30 
1  1  40 
1  150 
1  1  60 
1170 
1  1  SO 
11B1 
1182 
2000 
20  1 0 
2020 
2021 
2022 
2023 
2024 
2025 
2026 
2027 
2028 
2029 
2030 
203  1 
2032 
2033 
2034 
2035 
2036 
2037 
2040 
2060 
2061 
2062 
2063 
2064 
2080 
2100 
2120 
2140 
2219 
2220 
2221 

2240 
2260 
2280 
2290 


i  i  i  i  i  i  i  i  i  i  i  i  t  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  t  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 


i  i  i  i  i  i  i  i  i 


FROGRAM  :   THESIS: 


THIS  PROGRAM  IS  DESIGNED  TO  CONVERT  AND  NORMALIZE  THE 
THE  WAVEFORM  DATA  VALUE  AND  POINT  LOCATION,  TRANSFERED 
FROM  THE  NICOLET  3091  DIGITAL  OSC ILLASCOFE ,  TO  THE 
VOLTAGE  AND  TIME  VALUE.   THE  NICOLET  IS  COTROLED  BY 
THE  HP-3041  DATA  ACQUISITION  AND  CONTROL  UNIT. 
THE  FROGRAM  SQUARE  THE  VOLTAGE  VALUE  AND  THEN  TAKE 
THE  RUNNING  AVERAGE  FOR  FIFTY  SET  OF  READING.   TO 
INSURE  THE  DATA  IS  NOT  LOSSED  DURING  COMPUTER  ERROR, 
THE  DATA  IS  BEING  STORE  IN  THE  DISK  AFTER  EVERY  10 
SET  CALCULATION. 


I  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  r  i  I  I  i  I  I  I  I  i  i  i  i  i  I  I  I  i  i  i 


DIM  D ( 1 50 1 ) , T ( 1 50 1 ) ,3(1 50 1 ) 

DIM  VDATA* C20  3  ,VTIMEfC20]  ,  VCOUNT* CZOI! 

INTEGER  N, I 

INPUT  THE  FILE  NAME  TO  STORE  DATA 

DISP  "WHAT  IS  THE  FILE  NAME  FOR  VOLTAGE  DATA  ?" 

INPUT  VDATA-S 

DISP  "WHAT  IS  THE  FILE  NAME  FOR  TIME  DATA  ?" 

INPUT  VTIME-f 

DISP  "WHAT  IS  THE  FILE  NAME  OF  THE  COUNT  DATA  ?" 

INPUT  VCOUNT* 

DISP  "YOUR  INTEREST  TIME  INTERVAL  IS  FROM  CSTART]  TO  CENDT3  " 

DISP  "WHAT  IS  THE  STARTING  TIME  OF  MEASURING  IN  MICRO-SEC)  ?" 

INPUT  START 

DISP  "WHAT  IS  THE  ENDING  TIME  OF  MEASURING  (  IN  MICRO-SEC)  ?" 

INPUT  ENDT 

t 

!    CHECKING  THE  INTERESTED  TIME  INTERVAL  FROM  CRT  SCREEN 
i 

LET  N=l 
i 

!  IF  THIS  IS  NOT  THE  FIRST  SET  OF  DATA,  THE  COMPUTER  IS  GOING 

!  TO  READ  THE  DISK  FOR  PREVIOUS  DAT  INFORMATION,  AND  CARRY  ON 

!  THE  JOB 
i 

DISP  "IS  THIS  YOUR  FIRST  DATA  SET  OUT  OF  THE  50  DATA  SET?" 
DISP  "PRESS  0  IF  YES  ;  PRESS  1  IF  NO" 
INPUT  ANS 
IF  ANS=1  THEN  GOSUB  READ_DATA 

FOLLOWING  IS  THE  COMPUTER  RECEIVE  THE  DATA  FROM  3091  AND 
THEN  CALCULATE  THE  VOLTAGE  AND  TIME  VALUE. 

DIM  I*C5] ,D$C 20OO0: ,N* C40D 


CONTROL  10,3 
CONTROL  10,4 
CONTROL  10,2 


I! 
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r2"?2 

2  TOO 
2301 
2302 
2303 
2304 
2305 
2306 
2307 
2308 
2310 
2311 
2312 
2313 
2314 
2315 
2316 
2320 
2321 
2322 

4-  -J  .£.  J 

2324 

2340 
2360 
2380 
2390 
2400 
2420 
2440 
2441 
2443 
2445 
2446 
2447 
24  60 
2480 
2500 
2520 

2530 
2531 

254  0 
2560 
2580 
2600 
2620 
2621 
2640 
2660 
26SO 
2700 
2720 
2740 
2760 
2770 
2771 
2780 
2800 


OUTPUT  70?  : "0FN4" 

ABOVE  STATEMENT  IS  LET  THE  CRT  TO  "LIVE"  CONDITION 

FOLLOWING  ONE  IS  TO  LET  THE  CRT  IN  "STORE"  CONDITION 


OUTFUT  70?  ; "CLS5" 
WAIT  1000 

PRESS  BOTTOM  ONLY  USED  WITHOUT  THE  HP- 
DISP  "PRESS  3091  BUTTON" 


!421  UNIT 


LET  THE  DATA  TRASMIT  TO  THE  COMFUTER  THROUSHR  RS-232 

3UTFUT  70?  ; "CLS2" 

ABOVE  STATEMENT  LET  THE  WAVEFORM  INFORMATION  BE  TRANSFERED 
FROM  CRT  TO  COMFUTER  THROUGH  P.S-232  INTERFACE 


ENTER  10 
ENTER  10 
ENTER  10 


It 

Df 
N» 


HO=VAL  (NfC 16,20]) 
V1=(VAL  (Nf  C21  ,253)  -5!  MO 
H1=(VAL  CN*C31  ,353)  -5)  MO 


(VAL 
(VAL 


<NfC26.30J ) 
(NTC36,403) 


12) 

12) 


70? 

70? 

THE 

?0? 


OUTFUT 
OUTPUT 
'  LET 
OUTPUT 
J=l 
K  =  5 

PRINT  "N=". 
FOR  1=0  TO 
IF  I- START 


: "0PN5' 
; "0FN2' 
CRT  IN 
:  "CLS4' 


LIVE  MODE 


N 

ENDT 
OR  i: 


ENDT  THEN  2780 


!     NORMALIZE  THE  VOLTAGE  VALUE  AND  TIME  FOR  THE  INTEREST  INTERVAL 
■ 

D(I)=VAL  (Df CJ,K3) *V1 

T( I)=(I-HO) *H1 

D(I)=D(I)  "2 

IP  N=i  THEN  GOTO  2640 

S( I)  =  ( (N-l) »S( I)+D(  I)  ) /N 


N=l  THEN  S(I)=D(I) 
I=START  THEN  2760  ELSE  2680 
I-START+20  THEN  2760  ELSE  2700 
I=START+30  THEN  2760  ELSE  2720 
I=START+50  THEN  2760  ELSE  2740 
760  ELSE  2780 
":D(I):"VOLT    "j 
":D(I ) : "VOLT 


IF 

IF 

IP 

IF 

IF 

IF  I -ENDT  THEN 

DISP  T ( I ) ; "SEC 

PRINT  T< I ) ; "SEC 


SCI) 

"  :  S  ( I ) 


J-I+5 
l-=k+5 
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2920 

28  40 

2950 

2851 

2852 

2960 

2765 

2970 

2971 

2972 

2990 

2990 

3000 

3305 

3306 

3307 

3308 

3309 

3310 

3500 

3510 

3520 

3540 

3550 

3560 

3580 

3590 

3600 

3620 

3640 

3660 

3690 

3700 

3720 

374  0 

3760 

3  780 
3900 
3900 
3910 
3920 
3930 
3940 
3950 
4000 
4020 

4  0  4  0 
4060 
4061 
4080 
4  1  00 
4120 
4140 
4160 
4190 
4200 
4220 
4240 
4260 
4290 
4390 


NEXT  I 
N=N+1 


1    LET  THE  DATA  BEING  STORE  AFTER  EVERY  10  SET  OF  CALCULATION 

IF  N  MOD  10=1  THEM  2965  ELSE  2240 
IF  N=ll  THEN  2970  ELSE  2980 
CREATE  VDATA*, 200,8 

CREATE  VCOUNTr, 1 ,9 

CREATE  VTIMEJ-,200,8 
GOSUB  WRITE_DISK 
IF  N=51  THEN  3000  ELSE  2240 
END 


FOLLOWING  SUBROUTINE  IS  USED  7  0  WRITE  THE  DATA  TO  THE  DISK 


WRITE_DISK: 

FURGE  VDATAr' 
CREATE  VDATA* ,200.9 
ASSIGN*  1  TO   VDATA-* 
PURGE  VC0UM1 t 
CREATE  VCOUNTt, 1 ,8 
ASSIGNtt  2  TO   VCOUNTt 
FURGE  VTIMEJ- 
CREATE  VTINE* ,200,0 
ASSIGNtt  3  TO   VTIME* 
PRINT*  2  :  N 
FOR  I=START  TO  ENDT 


PRINT*  3 
PR  I NTH  1 
NEXT  I 

ASSIGNtt 
ASSIGNtt  : 
ASSIGNtt  ; 
RETURN 


T  (  I) 
5(1) 

TO  * 
TO  * 
TO   * 


FOLLOWING  SUBROUTINE  IS  USED  TO  READ  THE  DATA  FROM  PREVIOUS 
WRI TTEN  DISK 


READ_DAIA: 

ASSIGNtt  2  TO   VCOUNT* 

ASSIGNtt  1  TO   VDATA* 

READ*  2  ;  N 

PRINT  N 

FOR  I=START  TO  ENDT 

READ*  1  s  S(I) 

IF  I=START  THEN  4220  ELSE  4140 


IF  I=START+20  THEM  422' 

IF  I-START+30  THEN  422' 

IF  I=START+50  THEN  422' 
IF  [=ENDT  THEN  4220  ELG 
PRINT  S(I) 
NEXT  I 

ASSIGN*  1  TO  * 

ASSIGN*  2  TO  * 

RETURN 


ELSE  416' 
ELSE  4  19' 
ELSE  420' 
E  4  240 
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900 

BO  1 
B02 

803 
804 
805 
806 
807 
808 
809 
810 
81  1 
812 
817 
814 
815 
316 
880 
900 
910 
920 
930 
04,;, 

950 

960 

970 

980 

990 

991 

992 

993 

994 

996 

1 000 

1020 

1040 

1060 

1080 

1  1 00 

1120 

1  1  40 

1141 

1144 

1146 

1148 

1150 

1152 

1  154 

1  1  60 

1165 

1  1  70 

1180 

1 200 

1220 

1240 

1 260 

1280 

1290 

1291 


1  1  1  1  1  1  t  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  r  1  1  1  1  1  1  1  1  1  1  1  1  1 


PROGRAM  :  MAXIMUM 


1  1  1  1  1  1  1  1  1  1  > 11 


THIS  PROGRAM  13  DESIGNED  TO  READ  THE  SQUARED 
WAVEFORM  VOLTAGE  DATA  AND  TIME  DATA  STORED 
IN  THE  DISK  AND  THEN  CHOOSE  THE  LOCAL  MAXIMUM 
(WAVEFORM  ENSEMBLE).  THE  CHOOSED  MAXIMUMS  WILL 
BE  TAKEN  SUARE-ROOT  AND  STORE  THE  FINAL  RESULT 
ON  THE  DISK 


1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  i  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 


D I M  VDATA*  C  20  ]  ,  VT  I  ME*  C  20 1  .  VNUMBER*  [  20  I  .  VT  I  MEJ1AX*  C  20 1  ,  VMAX  I  MUMS*  C  20  ] 

DI3P  "WHAT  IS  THE  NAME  FOR  VOLTAGE  DATA  FILE  ?" 

INFUT  VDATA* 

DISR  "WHAT  IS  THE  NAME  FOR  TIME  DATA  FILE  "," 

INPUT  VTINE* 

DISP  "WHAT  IS  THE  NAME  FOR  FILE  TO  STORE  THE  NUMBER?  NO.  ?" 

INPUT  VNUMBER* 

DISP  "WHAT  IS  THE  NAME  FOR  FILE  TO  STORE  THE  TIME_MAX*  DATA  ?" 

INPUT  VTIME_MAX* 

DISR  "WHAT  IS  THE  NAME  FOR  FILE  TO  STORE  THE  MAXIMUMS*-?" 

INPUT  VMAX I  MUMS* 

DISP  "WHAT  IS  THE  STARTING  TIME  OF  YOUR  DATA  COLLECING  ?" 

INPUT  START 

DISP  "WHAT  IS  THE  ENDING  TIME  OF  THE  DATA  COLLECTING  ?" 

INPUT  ENDT 

1 

ASSIGNS  1  TO   VDATA* 
ASSIGNS  7  TO   VTIME* 
CREATE  VNUMBER* ,1,8 
ASSIGNS  4  TO   VNUMBER* 
CREATE  VTIME_MAX*-,  100,8 
ASSIGNS  5  TO   VTIME_MAX* 
CREATE  VMAX I  MUMS*, 100,8 
ASSIGNS  6  TO   VMAXIMUMS* 
1 

INTEGER  PRINTTER 

DISP  "WHAT  IS  THE  PRINTER  DEVICE  ADDRESS  NO.  ?" 
!   INPUT  THE  PRINTER  NO.   SCREEN  IS  1 
INPUT  PRINTTER 
PRINTER  IS  PRINTTER 
1 

DIM  S(1501) ,T(1501 ) ,TIM( 100) , MAXS ( 100) 
PRINT  USING  1170  ;  "TIME  (SEC. )", "VOLT" 
IMAGE  25X, 1 1A, 12X,4A 
INTEGER  I,  J 
FOR  I=START  TO  ENDT 
READS  3  ;  T ( I ) 
READS  1  ;  S(I) 
NEXT  I 
LET  J=0 
1 

1    IF  THE  FIRST  DATA  BEING  READ  IS  GREAT  THAN  THE  SECOND  ONE, 
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1292  '   THE  FIRST  ONE  IS  TREATED  AS  A  LOCAL  MAXIMUM. 

1293  ! 

1200  IF  S(START)>=  S(START+1)  THEN  1220  ELSE  1500 

1320      TIM (J )=T (START) 

1340      MAXS(J)=SQR  (S(START)) 

1360  PRINT*  5  ;  TIM (J) 

13S0  PRINTtt  6  ;  MAXSCJ) 

1400  PRINT  USING  1410  ;  " J=" ; J , TIM ( J ) ,MAXB ( J ) 

1410  IMAGE  12X , 3A , DDD , 7X , D. DDDDDD , 1 IX .D.DDDDDDDD 

1420  !  PRINT  TIM (J) 

1440  !  FRINT  MAXS(J) 

1460      J=J+1 

1 480  ! 

1500  FOR  I=START+1  TO  ENDT-1 

1520       IF  S(I).=  S(I-l)  AND  S < I )  >=  S  (  I ->- 1  )  THEN  1540  ELSE  1720 

1540  TIM(J)=T(I) 

1560  MAXS (J) =SOR  (S(I>) 

1580  ! 

1600  PRINT  USING  1610  ;  "  J  =  "  ,  J  ,  T  I  M  <  J  )  ,  MAXS  (  J  ) 

1610  IMAGE  12X, 3A , DDD. 7X,D. DDDDDD, 1 1 X , D.DDDDDDDD 

1620  !  PRINT  TIM (J) 

1640  !  PRINT  MAXS (J) 

1660  PRINTtt  5  ;  TIM (J) 

1680  PRINTtt  6  ;  MAXS(J) 

1700  J  =  J-t-l 

1720  NEXT  I 

1730  J=J-1 

1740  PRINTtt  4  :  J 

1760  ASSIGN**  4  TO   * 

1800  ASSIGN**  5  TO   * 

1820  ASSIGN**  6  TO   * 

1840  ASSIGN**  1  TO   * 

1880  ASSIGN**  3  TO   * 

1900  END 
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900 
901 

903 

904 

905 

906 

907 

90  B 

909 

910 

911 

912 

913 

914 

915 

916 

917 

919 

919 

1 000 

1 020 

1022 

1024 

1025 

1026 

1 028 

1040 

1041 

1042 

1043 

1060 

1080 

1  100 

1120 

1140 

1  160 

1130 

1200 

1220 

1240 

1241 

1242 

1243 

1244 

1245 

1  250 

1251 

1252 

1253 

1 1255 

1256 

1257 

1259 

1260 

1280 

1  300 

1320 

1340 

1360 


I  f  !  I  I  I  I  I 


PROGRAM 


PLOT" 


THIS  PROGRAM  IS  DESIGNED  TO  PLOT  THE  ENVELOF 
OF  WAVEFORM  BY  READING  DATA  FROM  DATA  DISK. 
AFTER  READING,  VOLTAGE  DATA  IS  CONVERTED  TO  LOG 
VALUE,  THEN  BE  PLOTTED.   AFTER  PLOTING,  A  LEAST- 
SQUARE  METHOD  IS  USED  TO  FIND  THE  BEST  FIT  STRAIGHT 
LINE  FOR  THE  LOG-ENSEMBLE,  PLOT  THE  LINE,  PRINT 
OUT  THE  SLOPE  OF  THE  LINE  AND  THE  CORRELATION 
COEFFICIENT. 


i  i  i  i  i  i  i  i  i  i  i  «  i 


i  i  i  t  i  i  i  i  t 


CLEAR 

GCLEAP 

INTEGER  PLOTER 

DISP  "WHAT  IS  YOUR  POLTTER  DEVICE  ADDRESS  NO.  ~>" 

!   INFUT  THE  PLOTTE  NO.   SCREEN  IS  1 

INPUT  PLOTER 

PLOTTER  IS  FLOTER 

FRAME 

THIS  PLOT  SIZE  IS  TO  FIT  THE  REQUIREMENT  FOR  THE  NFS  THESIS 

LOCATE  30,125,30,85 

CSIZE  3 

MOVE  50,20 

LABEL  "TIME  (MICRO-SEC)" 

PEN  UP 

MOVE  1 8 , 40 

DEG 

LDIR  90 

LABEL  "LN(  (  V0LTAGE'2>)  " 1 /2)  " 

PEN  UP 

DEFAULT  SET-UP  IS  TO  FIX  THE  PLOTING  COORDINATION  AS 
X-AXIS  IS  FROM  500  TO  670  (MICRO-SEC) ,  Y-AXIS  IS  FROM 
-5  TO  2 


DISP  "DO  YOU  WANT  TO  USE  DEFAULT  SET-UP 

INPUT  BULL 

IF  BULL=0  THEN  1253  ELSE  1260 

SCALE  500,670,-5,2 

LAXES  10,-5,500,-5,2,2 

XMIN=500 

GOTO  1710 
i 

DISP  "ENTER  THE  XMIN  OF  SCALE" 

INPUT  XMIN 

DISP  "ENTER  XMAX  OF  SCALE" 

INPUT  XMAX 

DISP  "ENTER  THE  YMIN  OF  SCALE" 

INPUT  YMIN 


YES=0' 
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1380  DIBP  "ENTER  THE  YMAX  OF  SCALE  " 

1400  INPUT  YMAX 

1420  SCALE  XMIN,XMAX,YMIN,YMAX 

1440  FXD  0,2 

1460  DISP  "ENTER  THE  X-TICKING  SPACE" 

14S0  INPUT  XT 

1500  DISP  "ENTER  THE  Y-TICKING  SPACE" 

1520  INFUT  YT 

1540  DISP  "ENTER  THE  X  INTERSECTION" 

1560  INPUT  XI 

1580  DISP  "ENTER  THE  Y  INTER  SECTION" 

1600  INPUT  YI 

1620  DISP  "ENTER  THE  X-MAJOR  COUNT" 

1640  INPUT  XMC 

1660  DISP  "ENTER  THE  Y-MAJOR  COUNT" 

1661  ! 

1680  INPUT  YMC 

1700  LAXES  XT, YT , X  I  , YI , XMC, YMC 

1710  GRAPH 

1720  PEN  1 

1721  ! 

1740  DIM  TIM(IOO) ,MAXS<100) 

1741  DIM  TIMEMAX*  [20  3  .MAXIMUM*  [20  3 

1742  DISP  "WHAT  IS  YOUR  FILE  NAME  FOR  TIMEMAX*  ?" 
1744  INPUT  TIMEMAX* 

1746  DISF  "WHAT  IS  YOUR  FILE  NAME  FOR  MAXIMUM*  7" 

1748  INPUT  MAXIMUM* 

1749  PRINT  "THE  DATA  FILES  ARE" , T IMEMAX* , " AND" , MAX IMUM* 

1750  DISP  "WHAT  IS  THE  HIGHEST  ORDER  OF  ARRAY  NUMBER'"" 

1751  '  THIS  GIVES  A  CHOOICE  TO  CHOOSE  THE  END  POINT  OF  PLOT 

1752  INPUT  JK 

1759  ! 

1760  ASSIGN**  5  TO   TIMEMAX* 
1780  ASSIGN*  6  TO   MAXIMUM* 

17B1  DISP  "PRESS  0  IF  THE  FIRST  LOCAL  MAXIMUM  IS  BEING  SAVED  ?" 

1782  INFUT  OK 

1790  FOR  1=0  TO  JK 

1800  READ*  5  :  TIM(I) 

1802  ! 

1803  !    CONVERT  THE  TIME  UNIT  FROM  SEC  TO  MICRO _SEC 

1804  ! 

1 S 1 0  T I M ( I ) =T I M ( I ) *  1 000000 

1820  READ*  6  ;  MAXS ( I ) 

1821  MAXS(I)=LOG  (MAXS ( I ) ) 

1822  !   CONVERT  THE  VOLTAGE  TO  NATURE  LOG  VALUE 

1823  NEXT  I 

1825  K=JK+1 

1826  IF  0K=0  THEN  GOTO  1835 
1828  FOR  1=1  TO  JK 

1831  TIM(I-1 )=TIM( I) 

1832  MAXSd-1  )  =MAXS(  I  ) 

1833  NEXT     I 

1834  K=JK 

1835  !  DISP  "WHICH  CHARACTER  YOU  PREFER  TO  USING  IN  PLOTTING  THE  DATA?1 
1837  !  INPUT  C* 

1839  MOVE  TIM(O) ,MAXS(0) 

1840  FOR  1=0  TO  K-l 
1920  DRAW  TIM( I ) ,MAXS( I ) 
1930  !  LABEL  C* 

1940  NEXT  I 
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1960 
1961 
1962 
1980 
2000 
2020 
2900 
29 1 0 
2920 
2922 
2924 
2926 
2928 
2930 
3000 
3040 
3060 
3060 
3 1  00 
3120 
3140 
3160 
3180 
3200 
3220 
3240 
3260 
3280 
3300 
3320 
3340 
3350 
3351 
3352 
3360 
3370 
3371 
3372 
3380 
3400 
3420 
3421 
3422 
3423 
3440 
3441 
3442 
3443 
3460 
3461 
3480 
3481 
3500 
3520 
3521 
3540 
3560 
3561 
3562 
3580 


PEN  UP 

6QSUB  LEASTSQR 
ASSIGN*  5  TO   * 
ASSIGN*  o  TO   * 
END 


FOLLOWING  SUBROUTINE  IS  DESIGNED  TO  APPLY  THE  LEAST 
SQUARE  METHOD  IN  THE  ABOVE  DATA  AND  THE  PLOT  THE 
BEST  FIT  STRAIGHT  LINE,  PRINT  OUT  THE  SLOPE  AND 
THE  CORRELATION  COEFFICIENT  ON  THE  PLOT. 
FORMULA  USED  PLEASE 

LEASTSQR: 

REAL  Y , XY , X 1 , X2 , TAVE , YAVE , VARX , VARY , COVXY 

LET  Y=0 

1         Y  IS  THE  SUMMATION  OF  /(I)  VALUE 
LET  XY=0 

!         XY  IS  THE  SUMMATION  OF  X < I ) *Y  (  I  ) 
LET  X1=0 

!         XI  IS  THE  SUMMATION  OF  X(I) 
LET  X2=0 

!  X2  IS  THE  SUMMATION  OF  XCI)  2 

FOR  1=0  TO  K-l 

Y  =  Y+MAXS  i I ) 

XY=XY+TIM< I) *MAXS(I ) 

X1=X1+TIM ( I) 

X2=X2+TIM( I)'  2 

NEXT  I 

B=(X2*Y-X1*XY> / <K*X2-X1  -2) 


A  IS  THE  SLOPE  OF  THE  FITTED  STRAIGHT  LIGHT 
A=(K*XY-X1*Y) / (K*X2-X1  2) 


B  IS  THE  INTERCEPT  TO  Y  AXIS 
PRINT  "THE  SLOPE  A=";A 
PRINT  "THE  Y-INTERCEPT  B=";B 
YAVE  = (MAXS  (0) +MAXS (K-l )  ) /2 

yave  is  the  average  value  for  the  log  voltage 
tave= ( t i m ( 0 ) +t i m ( k- 1 ) ) / 2 

tave  is  the  average  value  for  the  time  variable 
:ovxy=o 

covxy  is  the  covariance  of  time  and  voltage  <ln> 

VARX=0 

!        VARX  IS  THE  VARIANCE  OF  TIME  VARIABLE 


VARY=0 


VARY  IS  THE  VARIANCE  OF  VOLTAGE (LN)  VARIABLE 


DIM  YL(IOO) 
FOR  1=0  TO  K-l 
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3600 
3620 
3640 
3660 
3680 
3700 
3720 
3740 
3760 
3761 
3730 
3800 
5820 
5840 
5860 
5880 
3900 

:9io 

3920 

3940 
3960 
3980 
3990 
3991 
3994 
3996 
3998 
3999 
4000 
4002 
4004 
4006 
4020 
4021 
4040 
4060 
406 1 
4062 
4080 


C0VXY=C0VXY+ (TIM( I) -TAVE)* (MAX3 ( I ) -YAVE) 
i 

VARX=VARX+ (TIM ( I ) -TAVE) "2 

VARY=VARY+ (MAXS ( I) -YAVE)  ~2 

YL(I)=A*TIM(I)+B 

!        YL  IS  THE  VALUE  FITTED  ON  THE  STRAIGHT  LINE 

NEXT  I 

PEN  2 

MOVE  TIM<0) ,YL(0) 

FOR  1=0  TO  K-l 

DRAW  TIM ( I ) , YL ( I ) 

NEXT  I 

CORCOE=COVXY/SQR  (VARX*VARY) 

CORCOE  IS  THE  CORRELATION  COEFFICIENT 

FRINT  "CORRELATION  COEFFICIENT  IS  "; CORCOE 
PEN  I 
DEG 

LDIR  0 
i 

!      FOLLOWING  PRINT  OUT  THE  SLOPE  AND  COORELATION  COEFFICIENT 

!      ON  THE  PLOT 

i 

MOVE  XMIN+80, 1 
i 

!      LABLE  ALWAYS  STARTS  FROM  THE  XMIN+80  IN  X  COORDINATES 

LABEL  USING  4021  ;  "SLOPE  =",A 

IMAGE  7A,DDD.DDDD 

MOVE  XMIN+80, 2 

LABEL  USING  4061  ;  "CORRELATION  COEFF=" , CORCOE 

IMAGE  19A.DDD.DDDD 

PEN  UP 

RETURN 
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BOO 
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i  i  i  i  i  i  i  i  i  i  i  t  i  i  i  i  i  i  i  i  i  i  i 


PROGRAM 


i  i  i  i  i  i  i  i 


AVERAGE 


I  t  I  I  I  I  1  I  I  I  !  !  I 


THIS  PROGRAM  IS  DESIGNED  TO  AVERAGE 
THE  DATA  SETS  TO  GET  A  AVERAGED 
WAVEFORM  DATA. 

i  i  i  i  i  t  i  t  i  i  t  i  i  t  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  t  i  i  i  i  i  t  i  i  i  i  i  i  i  i  i 

DIM  VDATf C253 ,VDATA*C25D ,VDAT1*C25: 

DIM  V 1  ( 1 500 )  , V2 ( 1 500  >  , V3 ( 1 500  > 

DISP  "WHAT  IS  THE  ORIGINAL  DATA  FILE  ?   DAT**" 

INPUT  VDAT* 

DISP  "WHAT  IS  THE  NEW  DATA  FILE  ?   DATA**" 

INPUT  VDATA* 

DISP  "WHAT  IS  THE  AVERAGED  DATA  FILE  ?  DAT1**" 

INPUT  VDAT1* 

ASSIGN**  1  TO   VDAT* 

ASSIGN**  2  TO   VDATA* 

CREATE  VDAT 1$, 8, 1 10 

ASSIGN**  3  TO   VDAT1* 

N  IS  THE  NUMBER  HOW  MANY  DATA  SET  BEING  AVERGED 

DISP  "WHAT  IS   YOUR  AVERAGE  No.  N  ?" 
INPUT  N 

AVERAGE  ONLY  APPLY  TO  THE  SAME  ARRAY  NO. 

FOR  1=522  TO  621 

READ**  1  ;  VI  (I) 

NEXT  I 

FOR  1=522  TO  621 

READ**  2  ;  V2  (  I  ) 

NEXT  I 

FOR  1=522  TO  621 

RUNNING  AVERAGE  IS  ADOPTED 

V3( I )=( (N-l ) *V1 ( I >+V2(I) > /N 
FRINTt*  3  ;  V3  (  I  ) 

PRINT  OUT  FOLLOWING  VALUE  AS  A  CHECKING  PURPOSE 

IF  1=522  THEN  1500  ELSE  1400 

IF  1=553  THEN  1500  ELSE  1420 

IF  1=600  THEN  1500  ELSE  1440 

IF  1=621  THEN  1500  ELSE  1520 
PRINT  VI (I) ,V2(I) ,V3 ( I) 
NEXT  I 
END 
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APPENDIX  C 
DECAY  ENVELOPE  GRAPHS  FOR  AGGREGATE 

This  appendix  includes  the  graphs  of  the  decay  envelopes  for  reflection  from 
aggregate  water  interlace  within  the  glass  tank.    Figure  C.l  to  C.9  show  the  original 
decay  envelopes  obtained  from  different  transducer  positions.   The  decay  envelopes  are 
also  described  as  following  data  sets: 

data  set  42 

data  set  43 

data  set  44 

data  set  45 

data  set  46 

data  set  47 

data  set  48 

data  set  49 

data  set  50 
Figures  CIO  to  CIS  are  the  decay  envelopes  resulted  from  the  running  average 
method.    They  are  also  described  as  following  data  sets: 

data  set  DAT2 

data  set  DAT3 

data  set  DAT4 

data  set  DAT5 

data  set  DAT6 

data  set  DAT7 

data  set  DATS 

data  set  DAT9 

data  set  DAT  10 
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Figure  C.  1     Decav  envelope  for  reflection  from  aggregate/water  interface, 
transducer  was  30  cm  above  surface  (data  set  42). 
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Figure  C.2     Decay  envelope  For  reflection  from  aceregate/watcr  interface, 
transducer  was  30  cm  above  surface  (tlata  set  43). 
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Figure  C.3     Decay  envelope  for  reflection  from  aggregate/water  interface, 
transducer  was  30  cm  above  surface  (data  set  44). 
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Figure  C.4     Decay  envelope  for  reflection  from  aggregate/water  interface, 
transducer  was  30  cm  above  surface  (data  set  45). 
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Figure  C.5     Decay  envelope  for  reflection  from  aggregate/water  interface, 
transducer  was  30  cm  above  surface  (aata  set  46). 
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Figure  C.6     Decay  envelope  for  reflection  from  aggregate/water  interlace, 
transducer  was  30  cm  above  surface  (data  set  47). 


S2 


LN((<V0LTAGE-2>)~l/2) 


i 
en 


i 


n 

TO 

O 

1 

CD 

m 
r> 


500 


520 


540 


560 


580 


600 


620 


640 


660  - 


i 


i 
ro 


r\) 


CO 

o 

r~ 

o 

o 

33 

~a 

33 

m 

m 

r- 

ii 

> 

— i 

1 4 

i 

o 

o 

o 

n 

CO 

o 

en 

m 

~n 

Tl 

u 

CO 
CO 
CO 

o 


Figure  C.7     Decay  envelope  for  relleetion  from  aggregate/water  interface, 
transducer  was  30  cm  above  surface  (clata  set  48). 


83 


LN((<V0LTAGE"2>)"l/2) 


i 
en 


2: 

m 


o 

o 

I 

CO 

m 


500 


520 


540 


560 


580 


1 


I 
ro 


1 

i — • — r 


ro 


600 


620 


640  - 


660  - 


01 

D 

r~ 

O 

0 

;o 

"D 

:o 

m 

m 

r~ 

11 

> 

— 1 

» — • 

1 

0 

• 

:z 

0 

U) 

n 

•-J 

0 

^3 

rn 

T| 

00 

en 


Figure  C.8     Decay  envelope  for  reflection  from  aggregate/water  interface, 
transducer  was  30  cm  above  surface  (Uata'sct  49). 
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Figure  C.9     Dccav  envelope  for  reflection  from  aggregate  water  interface, 
transducer  was  30  cm  above  surface  (data  set  50). 
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Figure  CIO     Decay  envelope  resulted  from  the  average  of  data  set  42  and  43 
for  the  aggregate/water  surface  (data  set  DAI  2). 
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Figure  C.l  1     Decav  envelope  resulted  from  the  average  from  data  set 42  to  44 
lor  the  aggregate/water  surface  (data  set  DA  I  3). 
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Figure  C.12     Decay  envelope  resulted  from  the  average  from  data  set42  to  45 
for  the  aggregate/water  surface  (data  set  DA'14). 
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Figure  C.13     Decay  envelope  resulted  from  the  average  from  data  set 42  to  46 
for  the  aggregate/water  surface  (data  set  DAT5). 
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Figure  C.14     Decay  envelope  resulted  from  the  average  from  data  set  42  to  47 
for  the  aggregate/water  sin  lace  (data  set  DA'16). 
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'igure  C.15     Decay  envelope  resulted  from  the  average  from  data  set  42  to  48 
for  the  aggregate/water  surface  (data  set  DA  17). 
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Figure  C.16     Decay  envelope  resulted  from  the  average  from  data  set  42  to  49 
for  the  aggregate/water  surface  (data  set  DATS). 
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Figure  C.17     Decay  envelope  resulted  from  the  average  from  data  set  42  to  50 
for  the  aggregate/water  surface  (data  set  DA  19). 
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Figure  C.18     Decay  envelope  resulted  from  the  average  from  data  set  42  to  5 
for  the  aggregate/water  surface  (data  set  DAT  10). 
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APPENDIX  D 
DECAY  ENVELOPE  GRAPHS  FOR  FINE  SAND 

This  appendix  includes  the  graphs  of  the  decay  envelopes  for  reflection  from  fine- 
sand  water  interface  within  the  glass  tank.    Figure  D.l  to  D.10  shew  the  original  decay 
envelopes  obtained  from  different  transducer  positions.   The  decay  envelopes  are  also 
described  as  following  data  sets: 

data  set  55 

data  set  56 

data  set  57 

data  set  58 

data  set  59 

data  set  60 

data  set  61 

data  set  62 

data  set  63 

data  set  64 
Figures  D.ll  to  D.19  are  the  decay  envelopes  resulted  from  the  running  average 
method.   They  are  also  described  as  following  data  sets: 

data  set  DAT  12 

data  set  DAT  13 

data  set  DAT  14 

data  set  DAT  15 

data  set  DAT  16 

data  set  DAT  17 

data  set  DAT  18 

data  set  DAT  19 

data  set  DAT20 
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Figure  D.l     Decay  envelope  for  reflection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  55). 
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Figure  D.2     Decav  envelope  for  reflection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  56). 
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Figure  D.3     Decav  envelope  for  reflection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  57). 
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Figure  D.4     Decay  envelope  for  rcllcction  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  58). 
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Frigure  D.5     Decay  envelope  for  reflection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  59). 
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Figure  D.6     Decay  envelope  for  reflection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  60). 
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Figure  D.7     Decay  envelope  for  reflection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  61). 
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Figure  D.8     Dccav  envelope  for  reflection  from  fine-sand/water  interlace, 
transducer  was  30  cm  above  surface  (data  set  62). 
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Figure  D.9     Decay  envelope  for  reflection  from  fine-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  63). 
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Figure  D.10     Decay  envelope  for  reflection  from  line-sand/water  interface, 
transducer  was  30  cm  above  surface  (data  set  64). 
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Figure  0.11     Decay  envelope  resulted  from  the  average  of  data  set  55  and  56 
lor  the  fine-sand/water  surface  (data  set  DAT  12). 
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Figure  D.12     Decay  envelope  resulted  from  the  average  from  data  set  55  to  57 
for  the  fine-sand/water  surface  (data  set  DAI  13). 
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Figure  D.13     Decav  envelope  resulted  from  the  average  from  data  set  55  to  58 
for  tlie  fine-sand/ water  surface  (data  set  DAT  14). 
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Figure  D.14     Decav  envelope  resulted  from  the  average  from  data  set  55  to  59 
for  tlie  fine-sand/vvater  surface  (data  set  DAT15). 
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Figure  D.15     Decay  envelope  resulted  from  the  average  from  data  set  55  to  60 
for  the  fine-sand/water  surface  (data  sef  DA'l  16). 
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Figure  D.16     Decay  envelope  resulted  from  the  average  from  data  set  55  to  61 
lor  the  line-sand,  water  surface  (data  set  DAT  17). 
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Figure  D.17     Decay  envelope  resulted  from  the  average  from  data  set  55  to  62 
lor  tlie  fine-sand/water  surface  (data  set  DAI  18). 
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Figure  D.18     Decay  envelope  resulted  from  the  average  from  data  set  55  to  63 
for  tlie  fine-sand/water  surface  (data  set  DAI  19). 
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Figure  I). 19     Decay  envelope  icsultcd  from  the  average  from  data  set  55  to  64 
lor  the  aggregate/water  surface  (data  set  DA  I  20). 
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